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The Uniform Motion of a Sphere through a Viscous Liquid. 

By R. W. Bukgess. 



1. The problem of finding at any point in a viscous liquid the velocity 
due to the uniform rectilinear motion of a sphere was first attacked by Stokes 
in 1851 and an approximate solution obtained. From this was derived the 
well-known Stokes' law, that a force 6n(iWa is required to maintain a uniform 
velocity W of a sphere of radius a through a liquid whose coefficient of vis- 
cosity is p. This solution was proposed merely as a limiting form to which 
the distribution of velocities may be supposed to tend as W approaches the 
limit zero, but has nevertheless proved itself well in accord with experimental 
facts, at least as a close first approximation. The serious mathematical incom- 
pleteness of the solution has been recognized, however, and several attempts 
have been made to improve on it. Whitehead, in 1888, attempted to modify 
this solution by taking into account terms depending on squares of velocities, 
but found it impossible by his method to satisfy the boundary conditions at 
infinity. In 1893, Rayleigh showed that Stokes' solution would be accurate if 
certain additional forces were introduced, the magnitude of the required forces 
furnishing some indication of the error of the solution. In 1911, Oseen found 
a solution which is satisfactory at infinity ; it also gives a good approximation 
near the sphere if the velocity or the radius of the sphere is small. Shortly 
afterwards Lamb, using a different method, presented Oseen's results in a 
simpler form, but still subject to the same limitations. In 1913, F. Noether 
attempted to take into consideration the terms neglected by Stokes, but merely 
derived Whitehead's correction over again. A few months later Oseen took up 
the same problem and tried to apply the principle of his previous solution, but 
fell into the same error that he had carefully pointed out in Whitehead's work. 
In the following paper Oseen's results are obtained by a simpler method than 
that of either Oseen or Lamb. The method is generalized, so as to give a solu- 
tion that satisfies the boundary conditions at the sphere to any desired degree 
of approximation. A correction is determined which takes account of all the 
terms quadratic in velocities, and differs from the corrections of . Whitehead, 
Noether and Oseen in that it satisfies the boundary conditions at infinity. 
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2. Using cylindrical coordinates (r, z), let u, w denote the components of 

the velocity, p the pressure, d the density, v ( = ^ j the kinematic coefficient of 

viscosity, and V the potential of the impressed forces. The differential equa- 
tions for the motion of a viscous liquid, in the case when the Z-axis is an axis 
of symmetry, are* 



du , du , du dQ 
-df +U lk +W lte=-dr- +V 
dw dw , dw dQ „„ 

_ + „_ s _ + „, w= ^_ + , v%j 



K£> 



(i) 



where 



Eliminating Q from these equations, we obtain 

3 /3w , du , 3w\ d (dw dw dw 

■=- _ — hw, _ — k m; -= — I — ^-1 —= \-u-= \- id —= — 



\_d_ldw_ 
I dr\ dt 



dz\dt +U dr + W ^)-dr\-dt +U -dr- + W ^z 



This equation is simplified by making use of the equation of continuity f 



If we put 



du u dw _ 
dr r dz ~ ' 

u= ^-, w = — ~, (3) 

r dz r dr 



this relation is satisfied ; hence the velocities may be found by differentiation 
from this function 4 1 , known as Stokes' current function, and the problem is 
reduced to the determination of 4/. The curves of intersection of planes through 
thee-axis with the surface *= a constant, will be the lines of flow or stream 
lines ; that is, the lines indicating the direction of the resultant velocity at any 
point. This system of curves provides a convenient method of representing 
graphically the motion given by any function 4\ 

* Bassett, "Hydrodynamics," Vol. II, p. 244. f Bassett, "Hydrodynamics," Vol. I, pp. 8, 12. 
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Substituting in (2) the values of u and w given in (3), we have* 

where 

3r 2 r 3r 3s 2 ' 

The preceding equation applies to any case in which the motion of the 
liquid is symmetrical with respect to the s-azis. We desire to discuss in this 
paper the motion due to a sphere which has been moving for an infinite time 
with uniform velocity W along the s-axis, and whose center is at (0, Wt) at 
time t. Since the conditions at a point depend only on its position relative to 
the sphere, i^ is evidently a function of r and z — Wt only. If s'=s — Wt, we 
have 

34' _ d4> dz' d4> 34 1 _ 34 3s' „ 7 34' 

~dz ~ Wdz ~ W ' ~df~ dz' dt ~~~ W 

Equation (4) then becomes 

( u *. + {w - W )£—^-vd)d4,=0, (5) 

a form on which we shall base our further discussion. 

3. The problem we shall consider is the determination of the function 4- 
of r and z' which satisfies (5) and is such that the velocity of the liquid at the 
surface of the sphere is the same as the velocity of the sphere, and that the 
velocity of the liquid at an infinite distance from the sphere vanishes. This 
problem was first attacked by Stokesf in 1851 and solved, approximately, for 
the case in which the liquid is very viscous; that is, v large, and the veloci- 
ties sufficiently small so that the terms involving the products of velocities and 
accelerations are comparatively unimportant. On these assumptions the 
approximate form of (5) is 

vD-D^=0, (6) 

for which we may easily obtain Stokes' solution 



4 = i Wa*(^--— )sin 2 0, 



where a is the radius of the sphere, p= "vV+s 2 and 



r 



tan0=-. (7) 

z 



* Cf. Bassett, " Hydrodynamics," Vol. II, p. 259. 

t Cambridge Transactions, Vol. IX (1851) ; Sei. Papers, Vol. Ill, p. 55. 
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From this is derived Stokes' well-known formula F=6jt^Tfo for the force 
required to maintain the motion. This formula has proved so accurate for 
small velocities that the approximations by which it was obtained have not 
always been fully recognized. Stokes himself, however, proposed the solution 
(7) only as a limiting form to which the velocities may be supposed to tend as 
W approaches zero. 

In 1888, Whitehead * attempted to modify this solution to take account of 
the neglected quadratic terms. The attempt was successful in finding correction 
terms satisfactory as far as the differential equation (5) and the boundary con- 
ditions at the sphere were concerned, but these terms failed, as Whitehead 
himself noticed, to satisfy the condition that the velocity at an infinite distance 
vanishes. The correction terms for the velocities along and perpendicular to 
the radius vector were found to be 



3 Wa/a A'/fl* . a 



1= - - _( - -1 ) ( 7 + w + - + 2jsm 9 cos fl. 



(8) 



Vp Ap s ' 2p* 2p 
It is obvious that when p is infinite these terms do not vanish as required. 

In 1893, Eayleighf showed that Stokes' solution would be accurate if cer- 
tain additional forces were introduced; the ratio of the forces thus required 
to the viscous forces actually considered gives a means of estimating the accu- 
racy of the approximation. This ratio can easily be shown to be — £ , which 

v 

increases indefinitely with p ; Stokes' solution is therefore not valid at infinity. 
4. In 1911, Oseen| introduced a new approximate solution, valid at 

infinity. Consider the term (w—W) ^- in equation (5). Near the surface of 

OS 

the sphere it is true that w — W is very small, and the neglect of this term is 
justified, but at any considerable distance it is not; in fact, w — If is very 
nearly — W, not zero. Oseen shows that the ratio of the neglected terms to the 

terms giving the viscous forces is , which becomes infinite with z, no matter 

v 

how small we take W. On account of these considerations at a large distance 

from the sphere, we take as an approximate form of (5) 

(-W~-vd)d^=0. (9) 



* Quarterly Journal of Mathematics, Vol. XXIII (1888), pp. 143-152. 

t Philosophical Magazine, Ser. 4, Vol. XXXVI (1893), p. 354. 

X Arkiv for Matematik, Astronomi, ooh Fysik, Bd. 6 (1911), No. 29. 
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Oseen's method is not founded on a consideration of equation (5) or (9), nor 
does he use the current function ; but the above remarks embody his argument. 
He obtains a solution which, expressed as a current function, satisfies (9) and 
the boundary conditions at infinity, but satisfies the boundary conditions at the 

surface of the sphere only to the first power of the quantity , which is 

v 

assumed to be small. 

5. The essential points of this work of Oseen were restated by Lamb* 
shortly afterwards, with the use of a less intricate method of analysis. Lamb 
points out especially the significance of the result, showing that the equation 
(9), or rather the similar equations used by Oseen, represents more accurately 
than does Stokes' equation (6) the nature of the motion at a large distance 
from the sphere. Lamb shows, moreover, that Oseen's solution more nearly 
represents the wake which follows a moving object, in that the velocity preced- 
ing the sphere is markedly smaller than that following. 

Two years later F. Noether,f referring in his work to that of Oseen, deter- 
mined as a correction to Stokes' current function (7) 

^ 1= k^T Sin2 6 C0S fl^JW+ap+a 2 ). (10) 

Corresponding to this, the corrections for the velocities are 

R >= - i6 — p *< cos 0> ( V A 2 + p + 7> 

@1= 16 -V Sm ° C ° S 6 V- p A 2+ 2p + 2? + 7 

which are exactly (8), as found by Whitehead. 

In the same year Oseen $ again published an article on the subject criti- 
cizing certain portions of Noether's work, and giving as the value of Stokes' 
current function 



3 Wa 1 W 

+ |I^cos0 (2p 2 +ap+a 2 ) j-^fsitfB. (11) 

But this function yields values of R and © which fail to vanish at infinity, 



* Philosophical Magazine, Ser. 6, Vol. XXI (1911), pp. 112-121. 

t Zeitsohrift fiir Mathematik und Physik, Bd. 62- 

% Arkiv for Matematik, Astronomi, och Fysik, Bd. 9, No. 10. 
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exactly as Whitehead's resolution did, and therefore, as Oseen pointed out in 
his earlier paper, is not admissible. 

6. We shall now consider in detail equation (9), 



(— W-^-vD)D^=0. 



(9) 



As these operations are commutative, a solution will be given by 

^=i'+V, (12) 

where 

D# = (13) 

and 

(vD+W-~)v=0. (14) 

We shall obtain the solution of these separately. Equation (13) has been 
treated very fully by Sampson,* but for convenience of reference we shall here 
give in detail the solution so far as is necessary for our purposes. Transform- 
ing to polar coordinates 

p=-vV+s' 3 , = cos- 1 -7-=- 2 =sin~ 1 -j==% , 
our operation D becomes 

8 2 sin 6 d I 1 d 
9p 2 + ,p 2 ddKsmOdd 

We desire, then, the solution of the equation 



A particular solution is 



a 2 4 sine d/ i a*\ 



. sm0 3P„ t-iow 

+«=^¥» (13A) 



where P„ is the zonal harmonic of degree n, and therefore satisfies the differ- 
ential equation 

n(n + l)P n sin 6+ ~(sin6 ~^)=0. (15) 

For substituting in (13) the given value of -^ n , we have for the left-hand side 



n(n+l)c n „dP n c n sm$ d 
sin v -^- -\ 



p »+2 "~ d$ T- ?n+2 d Q 



1 d / . a dP, 
— • — * ^ sin 



L sin 6d6\ 36 



* Philosophical Transactions, Vol. CLXXXII (1891), p. 449. 
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that is, using (15), 

C -^[n(n+l) d -^-^\n(n+l)P n \], 
which is identically zero. It is well known that 

P M =^icos» ft- ifc- 1 ! nno *^* , 



2(2n — 1) 



n(n-l) (n-2) (n-3) -1 

P+ 2 -4(2n-l) (2n-3) * ••■.)■ (16) 



Hence, by (14), 



^ = c m sin 2 e r^„»-^ (n-l) (n-2) 



- cos"- 1 ©- 

p" L 



2(2n— 1) 



cos"- 3 0+ 



(w— 1) (n— 2) (n—2r) 



+ (~l) f »H°/n^ v i\ - ' )o' V '"o t^ eos-»+ 1 fl+ • • • -1 (13B) 
2 r \r(2n — 1) .... (2n— 2r+l) J 

for all positive integral values of n. For n—0, we have ^ a function of only; 
then 

1 Wo I A 

Placing w successively 1, 2, 3, in (13B), we have 

sin 2 



4>i = Ci- 



P 



, sin 2 6 cos 6 

^ 2 = c 2 - 2 , 



^3 = Cj 



P 

sin 2 fl (5 cos 2 0— 1) 
P 3 



, sin 2 cos 0(7 cos 2 0—3) 
4x=c t £ >- 



(17) 



7. We desire also the solution of equation (14), which is, if 2k= — , 

v 



Letting fy' = e~ hz 'ty(r, z'), we have at once 

/a 2 i d , a 2 . \ _ n 
Ur 2 " 7a7 + aP - *7*-°' 

or, in polar coordinates, 

a> , sinoa/ i a*\ 
3p 2+ p 2 ae\iSl"W/ * 



(14A) 



(18) 
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Assume $=B®, where B and © are respectively functions of p and only. 
^B , o sin0 d ( 1 dS\ wn . . 

R dp' k p - 80\sin0 dd )- n ( n + 1 )> sa ^ 
We have, then, to determine R and ©, the two equations : 

d ^-m- n -^+^-R=0, (19) 

dp p 

.(. + i)e+rf.»K^,^)=o. (20) 

In this form, for w any positive integer, i? is a solution in finite form of 
Eiccati's equation ; that is, 

R n=p «+i(l^y +1 (Ae*>+Ber»), (21) 

dP 

and © is sin -*£ • We have, therefore, 

n=<a 3D 

<$" = 2 e- k > cose amO^R n . (22) 

»=o <W 

We have, then, as a general function satisfying the differential equation (9) 

^ =e ~^l+cosS)V C + J) cog Q + bi gin 2 gA + 1\ 

+ & 2 sin 2 dcos^ 2 + — +4) + ....] 

_ „ , Msin 2 L lysin 2 cos , Psin 2 2fl ,„. 

+Lcos 0-| 1 s 1 5 — (5 cos'' 0—1) + (Jo) 

9 9 9 

8. To adapt this solution to the problem of the sphere moving with uni- 
form velocity W, we must have, if 



1 di> . 
R = -t — : — 7i-^7r is the radial velocity 
p 2 sin dd 

1 d4> 

and @ — ; — *-^— is the velocity perpendicular to this, 

p sin dp 

(1) B and © always finite, 

(2) B = and 0=0 at infinity, 

and (3) B = W cos 0, ©= — W sin on the surface of the sphere, 

i. e., when p = a. 



(24) 



(25) 
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The terms we have chosen are all such that R and © derived from them vanish 
at an infinite distance, and, if we take C-\-D cos 6=b (l — cos 6), are everywhere 

finite, since -^- and -~— are both divisible by sm 6. The last condition [25 (3) ] 
ou op 

remains to be satisfied by a proper combination of the functions and determi- 
nation of the constants. 

Oseen's solution is equivalent to taking the terms in (23) whose coefficients 
are b , L, and M ; that is to say, 

^=& o e- w < 1+cos «(l-cos0)+Lcos0+^^. (26) 

Using the conditions [25(3)] at the surface of the sphere to determine the 
constants, we have 



w L , 2i¥cos0 , b e- ka < 1+co ° e) M , . M fl . , 

Wcos 0= — -j- + -3 — + -s \1 + ka (1— cos d)\, 

. W sin e= ^i^ + M e-Mi+cos*) gin 0> 



(27) 



Equating the constant terms in each of these two equations and the coefficients 
of cos 6 on each side of the first, we have 

L -h TV- 2M 2fe o fe w- M b ° k 

a a a a a a 

It is to be noticed that we have neglected all squares and higher powers of the 
quantity ha, which must therefore be assumed to be small. Hence we derive 

T . Zav ,. Wa 3 

L = b Q = 2~, M = 4~- 

Substituting these values in (26), we have 

^ = _ §y [e -^a + ».«) ( i_ eo8 g) +cos 6]- WaS * in * e . (28) 

If we transform this value of ^ into cylindrical coordinates, find u and w from 
equations (3) , and resolve u in the directions of the x and y axes, we get expres- 
sions for the velocities which, except for notation, are exactly those derived 
by Oseen. As explained in section 4, this solution differs from Stokes' solu- 
tion (7) fundamentally in that it satisfies the differential equation (9), which 
for large values of p is a close approximation to the complete equation (5). 
We can see its significance by throwing it into several different forms. In the 
first place, since a constant makes no difference in a function used solely in the 
12 
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form of its derivatives, we may write it 

^ _ 4 L — Ml + cos0) pJ 
Expanding the exponential term, this is 

^^g_!) stf0 - 3 -^^in^(l+co S fl) + (28B) 

The first term of this expansion is Stokes' solution; the following terms con- 
tain powers of *- and are therefore less important than the first when p is 

(A. 

small (i. e., near the sphere) and W is small (i. e., for slow motion). As 
Oseen says, then, in the neighborhood of the sphere the solution approaches 
that of Stokes as a limit when W approaches zero. It is important to notice, 
however, that the character of the s'olution is essentially different from that of 
Stokes. If we plot the stream-lines for Stokes' solution, we notice that we 
have symmetry with regard to the plane z' =0 ; for Oseen's, however, the stream- 
lines are more crowded following the sphere. If, for instance, we set the first 
two terms of (28) equal to a constant; that is, e - k " a+COB6) (l— cos 0) +cos 6=c, 
the curves thus resulting, more conveniently expressed in the form 

. 1— cos 

_ l0g F=^sT 

p &(l + cos0) ' 
are asymptotic to the lines 0=cos _1 c before the sphere and run off in a para- 
bolic form in the wake of the sphere. The simplest way of seeing the differ- 
ence between the conditions before and after the sphere, is to notice that the 
velocities will each have an exponential factor — Zcp(l + cos 0) ; when is in the 
first quadrant, this term will be smaller than when is in the second quadrant, 
for the same values of p. Since this corresponds with our observations of fluid 
motion, the solution derived by Oseen appears worthy of further consideration 
and improvement. 

Expanding the boundary conditions as given in (27), we have 

^ . L 2Mcos0 ^T , ... ., k 2 a? ... „,„ 

Wcos6= j- + - 3 — - + -!U-Jfca(l+cos0) + -— (l+cos0) s 

Ot & Of L I ^ 

x{l+ha(l— cos0) j, 



\2 



M sin bji sin 



-W sin 0= ; + 



a" a 



2_2 



1— &a(l+cos 0) + -r^- (1+cos 0) 



(27A) 



From these equations it is clear that on the surface of the sphere the velocity 
components found in (28) differ from the velocity components of the sphere 
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by terms whose ratio to the terms W cos and — W sin 0, is ha. We proceed 
to improve this approximation by using two more terms from (20), and by this 
means making the above-mentioned ratio equal to h 2 a 2 . 
9. Let 



4=sin 2 0^ o ^— g- +& 1 (fc+-je 

- 1+COS0 \ p 



-io(i + co S v + M + Ncos6 



■■ sin 2 



P P 



(29) 



+ 



+ k(=^ + li^(, + i)l( 1+C os0r + ....](29A) 

L | w+1 [w \ p/J J 

= sin 2 0[5 o +£i(1 + cob0) +# 2 (l + cos0) 2 + .... +£ ra (l + cos0)" + . . . . ].(29B) 
The boundary conditions (21) and (22) at the surface of the sphere give 

1 8^ 



R. 



= — 



p 2 sin0 30 

l 34 



■-WoosB, 

= — W sin 0, when p=a. 



p sin dp 
As the first of these is to be true, for all values of 0, we may integrate, giving 

4= -Wa? sin 2 0, when p=a. We have, therefore, 



±Wa? sin 2 0=sin 2 [£„+#! (1+cos 0) + . . . . ], 
Wa sin 2 0=sin 2 0^ + f^J (i + cos 0) + ... .]. 



Dividing by sin 2 and equatiug like powers of (1 + cos 0), we obtain 



-Wa % =B, 



0=B n ] t 



•] - 

Jp=o 



Op Jp=a 
3p Jp=a 



(30) 



Using the values of B , B lf etc., given above in equation (29A), we find that 

37? 37? 

we can satisfy the four equations for B ,B t , -~ and -^ , neglecting squares 

and higher powers of ha, with the following values for the constants : 
b -- 



h=- 



\wa*{l-\ka), 



.. Wa'L 3 7 , 



N: 



-Wa A -ha. 
o 



(31) 
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To satisfy the boundary conditions (30) exactly, we should have B n and 

37? 

-^- 2 =0 when p = a; it is obvious that, while the given values of the constants 

do not effect this, the terms which should but do not vanish all contain at least 
the square of ha as a factor, and our error in approximately satisfying the 
boundary conditions at the sphere is proportional to the square instead of the 
first power of the small quantity -ka. We see, therefore, that while solution 
(28) holds as valid in the limit when ha approaches zero, the solution just 
found is a better approximation when ha is sufficiently small so that its square, 
but not its first power, may be neglected. If a closer approximation is desired, 
two more terms may be added, namely: 

7 »»nj.«,„* • <>n Jti , 3fc , 3\ , P sin 2 0(5 cos 2 0—1) ; ocy , 

V p fl f 

and we shall then be able to satisfy exactly the boundary conditions B 2 =0, 

37? 

-p- 2 - =0 when p=a, making our error proportional to h 2 a s W. Similarly, we can 

add any number of pairs of functions and make our error proportional to as 

high a power of ha as we desire. On account also of the jj^in the denominator 

of the n-th. term of (29), this error is more and more negligible as we go on. 

We have shown, then, how to obtain as accurate a solution of equation (9) as 

'i 
may be desired. 

10. In this differential equation (9) which we have just discussed, we have 

neglected certain terms which were included in our fundamental equation (5). 

We shall now include these and proceed to obtain a solution of the complete 

equation (5) accurate as far as the terms involving h?a 2 or /t 2 p 2 as a factor. 

When transformed to polar coordinates, equation (5) becomes 





1 

sin 


(ty 
0V30 


3 
3 P 


3^ 


d\, 


2 


; 0l COS V 

3 \T3 2 4 

30/ JL 3 P 2 


sin 34- \ 
p 30/ 

f 3 2 sin 3 / 1 
"l3p 2 ' p 2 30V sin 
sin 3 / 3^ \"| 
+ p 2 30\sin030/J - ' 




P 2 


3p 30/ ' p 2 
Wcos0^ — 


sin 2 

sin 
P 


3' 

'0 30, 

0. ( 



(5A) 

Let us assume that the desired approximate solution of this equation is of the 
form 

4=sin 2 0[A+Acos0], (33) 

where A and A x are functions of p, and A x contains ha as a factor. This 
assumption appears appropriate in view of the fact that, to the desired degree 
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of approximation, each of the solutions (7), (28) and (29) already found is 
of this type. Substituting this value of 4> i n tne differential equation (5A), 
omitting all terms which contain k 2 a? as a factor, we have, after some simplifi- 
cation, 



(f-mw-fyM-Th-* 



raM, 



4 d 2 A» 8_ 3A Q 8A 

"T .a 



3p 4 p 2 3p 2 ' p s 3p p 4 

. (&A X 12 3% , 24 3^A J n /wv 

+(-37— 7^+7¥/ C0 J (34) 

In order to satisfy this relation, we must in the first place have 

d'Ap ±^A 8 L dAo_S_A, .„_> 

3p 4 p 2 3p 2 + p 3 3p p 4 ' v ; 

That is to say, any value assumed for A must be such that A sin 2 $ satisfies 
D 2 '4'=0, since (35) results from this part of (5A), which was used by Stokes. 
The complete solution of this equation can easily be shown to be 

A = — +a 1 p + a 2 p 2 +a 8 p 4 - 
In order to satisfy (34), we must also have 

Hence we may easily determine the value of A x corresponding to any value 
of A . 

11. The first obvious assumption is to take A from Stokes' solution (7), 

'3p aN 



Then 



a p. 



&A 2A 6b 



(37) 



3p 2 p 2 ap ' 

and (34) becomes 

,^- 1V <^ +2 ^ = ^-^,). (38, 

Assume a solution, 

Al =d l9 * + d 2 p+d s +^ + ^, (39) 

in which the third and fifth terms are parts of the complementary function, 
the fifth term being essentially ^ 2 of equation (17). We could now easily 
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obtain the value of 4; our work is simplified, however, by first determining 
b . As our solution must satisfy the boundary conditions at the surface of 
the sphere, we have, as in (25), 



d±l 



w cos 0, 



1 d+i 

sin dp Jp=a~ 



a 2 sin 6 30 j P = a 
Since 4'=sin 2 0(4+4 cos 0), 

Wo? cos 0=24, cos 0+4 [—1+3 cos 2 6], 



W sin 0. 



(40) 



„ 7 dA 34 „ 
Wa= -^ + -^cos 0. 



We must therefore have 

Wa* 



dp dp 

=4] , 



Wa= 3 -f\ , 

Op Jp=a 



li] =0, 



dp Jp=o 



(41) 



From these conditions on A , we find b = jWa 2 . Substituting this value for 
b , (38) becomes 

„4 ^ 



10 ,3 2 4 , _ . 34 
P""87"~ 12 P^ +24 P 



¥ 



3p 



9TF 2 a 
4vp 



[— 2p 8 +3ap 2 — a 8 ]. 



Hence we easily obtain 



4=- 

Using (41), we have 



3 W*a 



32 y 



[2p 4 — 3ap s +d s p 2 a 2 — a s p+d 5 ] 



d s —d R =l. 



(42) 

(43) 
(44) 



. 3 W*a L a\ 2 /n „ , , „. 

•• ^ = -32~ ( r( 1 -p/ (2 P +a P +a) - 

The solution given by this assumption of the value of A is therefore 



*=!*.*> fa-?)-*" 



8 »1— r ) (2p 2 +ap+a 2 ) cos 



(45) 



which is exactly the correction found by Noether, as stated in (10) ; the veloci- 
ties derived from this are given in (8), and were first found by Whitehead. 
The objection which Whitehead himself raised against this approximation still 
stands, however ; that is, the velocities do not vanish at infinity as required by 
the conditions of the problem. 

12. In order to avoid this difficulty at infinity, let us take A from the 
solution of an equation which is, as Stokes' is not, a good approximation at a 
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large distance from the sphere. The simplest one of this type is (28). We 
shall consider this as our basis and determine a correction, accurate as far as 
terms involving k 2 a 2 or k 2 p 2 as a factor, which shall take account of the quad- 
ratic terms. Accordingly we take 

Determining b at once from equations (41), we have 



Wa 2 

2 


= 2& — \kaWa 2 , 


Wa= 


■■ — - — -kaWai 
a 4 


6o= 


ra+i*.). 



that is 



Substituting the value of A from (46) in (36), and neglecting, as we are doing 
throughout, all terms involving the square of ka, we again obtain (42). The 
approximate solution ij/=sin 2 B{A -\-A 1 cos 6) derived by this method is then 

*= \ W sin 2 e[(l+ \ka){z 9 a- |)-^V 

— |fta -cob (l— -Y (2p 2 +ap + a 2 )], (47) 

which is in substance exactly equation (11) , as derived by Oseen. In this form, 
however, the solution is still subject to the objection that it gives velocities 
which do not vanish at infinity. 

Solution (47) is therefore unsatisfactory; but if we notice that the cor- 
rection taking into account the quadratic terms is based on solution (28), and 

must therefore be added to it to give the required result, this difficulty disap- 

3 
pears. The term — ^ Wkap 2 sin cos of (47) is already included in the expo- 
nential term of (28) ; combining with (28) the terms of (47) not already 
included, we obtain 



^--3^(l+|jfca)sin 2 



g— kfi(l+cos6 



1 + COS0 

Wi *(i _■_ 3 i. \ sin2 e _i ' SWka • 2 a J* 2^« 3 « 4 \ t±a\ 

-r\}+l ha )-y- + ~iq- sm e cos e(3ap-a 2 + j - yj. (48) 
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Hence, as usual, 



p 2 sin0 80 " 2p 
p sin dp 



(l+.foj +__^_ 1+3 cos 2 0)(^ 3ap _a 2 +-_ j,( 49) 



sm0 8p 4 \ 4 / 



Wa 



A3, \sin0 3TfA;a . , J\ a' 2aH ,- m 
(l+j^j-3 jg- sin cos 0|_3a- -p + yj- (50) 



4 \ ' 4 / p 8 16p L p 2 ' p 3 

Using these velocities, we find with Oseen that Stokes' law should be 

F = 67t(iWal 1+ — g — ). 

As is obvious, these velocities, (49) and (50), both vanish when p is infinite; 
when p=a, they are, except for terms involving the squares or higher powers 
of ka, equal to W cos and — W sin 0, respectively ; and since, omitting squares 
and higher powers of ka and kp, (48) is the same as (47), to that approxima- 
tion they satisfy the complete differential equation (5). The solution, then, is 
superior to any previously given. The difficulty first raised by Whitehead, 
that a certain term gives velocities which do not vanish at infinity, is met by 
showing that this term is included in an exponential expression. 

Cornell University, May, 1915. 



